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Abstract. Low energy hadron phenomenology involving the (u,d,s) quarks is often ap- 
proached through effective multi-quark Lagrangians with the symmetries of QCD. A very 
successful approach consists in taking the four-quark Nambu-Jona-Lasinio Lagrangian with 
the chiral Ul(3) x Ur(3) symmetry in the massless limit, combined with the ?7a(1) breaking 
six-quark flavour determinant interaction of 't Hooft. We review the present status and some 
very recent developments related to the functional integration over the cubic term in auxi- 
liary mesonic variables that one introduces to bosonizc the system. Various approaches for 
handling this functional, which cannot be integrated exactly, are discussed: the stationary 
phase approximation, the pcrturbative expansion, the loop expansion, their interrelation 
and importance for the evaluation of the effective action. The intricate group structure 
rules out the method of Airy's integral. The problem of the instability of the vacuum is 
stated and a solution given by including eight-quark interactions. 
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1 Introduction 

Non-perturbative methods to deal with the low energy regime of QCD involve lattice calcula- 
tions, QCD sum rules pP, effective field theories, which are renormalizable in the modern sense 
and gave birth to Chiral Perturbation Theory dHHj, large N c (number of colours) expansions [IJG3 
HHEj, and effective models with solid symmetry content jHj. The object of our studies is placed in 
the latter category, consisting in a model Lagrangian which combines the four-quark interactions 
of Nambu and Jona-Lasinio (NJL) jHj, extended to the Ul{2>) x Ur (3) chiral symmetry of massless 
QCD ^3^2E1> and the six-quark interaction of 't Hooft which breaks the unwanted J7a(1) 
symmetry of QCD. We call it the NJLH Lagrangian. The 2Nf (Nf denoting the number of 
flavours) multi-quark interactions of 't Hooft were derived after eliminating the gluonic degrees 
of freedom of QCD in the semiclassical instanton approximation, assuming dominance of the 
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zero modes in the scattering process of quarks and anti-quarks on the same instanton. Recent 
reviews providing further evidence on multi-quark vertices can be found in |14j . 

In the last decade some new developments brought further relevance to the NJLH model, 
as it may help to shed some light on the possible existence of an hierarchy of multi-quark 
interactions, with prevalence of lower ones. There is some conflict on two fronts: (i) Going 
beyond the zero mode approximation, an infinite number of multi-quark interactions, starting 
from the four-quark ones, has been found in the instanton-gas model, all of them being of the 
same importance [T^]- (p) On the other hand, accurate lattice measurements for the realistic 
QCD vacuum show a hierarchy between the gluon field correlators with a dominance of the lowest 
one jnj- A similar hierarchy of the multi-quark interactions can be triggered after averaging 
over gluon fields. If this is true, there is an apparent contradiction with the instanton-gas model 
(see also |17j). This point has been stressed in |18j . 

The hierarchy problem of multi-quark interactions can be addressed on pure phenomeno- 
logical grounds. For that we suggest to simplify the task considering only the four- and six- 
quark interactions of the NJLH model. Once one knows the Lagrangian the obvious question 
arises: does the system possess a stable vacuum state and does this state correspond to our 
phenomenological expectations? If hierarchy takes place this question is pertinent for the leading 
four-quark interaction, because in this case the effective quark Lagrangian can be studied step 
by step in the hierarchy with the assumption that four-quark vertices are the most important 
ones. In the opposite case we must study the system as a whole to answer the question. For the 
best known and simplest example, to which we dedicate most of our attention here, the solution 
may be found analytically. As a result one can obtain definite answers to the above questions 
with a convincing indication in favour of a hierarchy for the considered example. 

Pioneering work on the original NJLH Lagrangian has been done in |19j . using Bethe-Salpeter 
techniques to evaluate quark - anti-quark scattering in one-loop approximation and obtain the 
meson masses from the poles of the scattering matrix, and in |2U| . where the quark Lagrangian 
has been bosonized using functional integral methods and the mesonic spectrum calculated at the 
leading order of the stationary phase approach (SPA) for one critical point. A good description of 
the pseudoscalar nonet, especially the 77 and rj masses and mixing has been achieved. Within the 
respective approximations, the effective potentials for both methods coincide The model has 
subsequently been widely and successfully explored at the mean field level, see e.g. [22 } l23l 124"] . 

In the present work we present a detailed picture and mathematically careful study of the 
functional integration over the auxiliary bosonic fields of cubic order, inherited from the 't Hooft 
six-quark interaction. They cannot be integrated out exactly. We came across one startling fact 
and several interesting properties, which we hope help in the understanding of the approxima- 
tions used and in the identification of the shortcomings and potentialities of multi-quark models. 
Before turning to the calculations, let us summarize what we have learned from them |25l I26j : 

(1) We show that the stationary phase equations have several roots (critical points), of which 
one is regular and the others singular. A rigorous SPA treatment, taking into account all critical 
points, leads to an unstable vacuum for the theory. This is the startling point, since this fatal 
flaw of the model is in crass contradiction with its phenomenological success. 

(2) The results in [2U] are obtained if one considers only the regular critical point of the SPA, 
which of course is accompanied with a different asymptotic behaviour. 

(3) This result can be understood in the perturbative sense, where the cubic interaction is 
taken as a perturbation around the stable NJL vacuum. The situation is in many respects 
analogous to the problem of a harmonic oscillator perturbed by an x s term. This system has no 
ground state, but perturbation theory around a local minimum does not know this. 

(4) By performing a resummation of the perturbative series in powers of the 't Hooft inter- 
action one gets a loop expansion [2Zj , of which the leading term contains all tree graphs present 
in all powers of the perturbative series. This corresponds exactly to the result of [2U] . 
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(5) The loop expansion can be classified in terms which contribute to the measure (all odd 
numbers of loops) and to the phase, relevant for the effective action (all even number of loops). 

(6) The global instability of the vacuum can be removed within the multi-quark interaction 
picture: the addition of a chiral invariant OZI-violating eight-quark interaction to the 
Lagrangian renders it globally stable, for certain values of the coupling strengths of the several 
interactions. 

(7) A comparison between the effective potentials for the NJLH and for the eight-quark 
enlarged Lagrangian reveals that the global stabilization happens around the local minimum 
which appears if only the regular critical point is considerd in the NJLH Lagrangian. This is 
probably at the heart of its success. 

Points (6) and (7) will not be discussed in this contribution. In fact, at the time the conference 
took place, they had not been obtained yet. We added them to the list, as they represent a kind 
of "happy-end" for the model. We refer to [201 f° r a detailed presentation of this issue. 

The paper is organized as follows. In Section 2 we write out the NJLH Lagrangian and use 
the functional integral representation to express it in terms of bosonic fields. In Section 3 we 
show that the chiral symmetry group imposes constraints which are only compatible either with 
the perturbative approach, or the expansion in a parameter that multiplies the total Lagrangian 
density (the loop expansion). Otherwise, as the consistent stationary phase treatment shows, 
the model is unstable. In Section 4.1 we discuss the perturbative treatment; the integration over 
auxiliary variables leads to a special problem with 6(0) -singularities, to which we give a physical 
meaning through Feynman diagrams in 4.2. The loop expansion is considered in Section 5. We 
obtain in closed form the two-loop contributions to the functional integral Z and give arguments 
to justify this result. Since we have forestalled the conclusions in the Introduction, we shall omit 
concluding remarks. 



2 The Lagrangian and bosonization 

On lines suggested by multicolour chromodynamics it can be argued |29| that the XJpSX) anomaly 
vanishes in the large N c limit, so that mesons come degenerate in mass nonets. Hence the leading 
order (in N c counting) mesonic Lagrangian and the corresponding underlying quark Lagrangian 
must inherit the Ul{3) x Ur(3) chiral symmetry of massless QCD. In accordance with these 
expectations the Ul(3) x Ur(3) symmetric NJL interactions, 

G 

£njl = [{qKqf + {qiisKq) 2 } , (l) 

can be used to specify the corresponding local part of the effective quark Lagrangian in channels 
with quantum numbers J p = + , 0". The Gell-Mann matrices and the singlet Ao = y| 1 act in 
flavour space and are denoted by A n , a = 0, 1, . . . , 8, obeying the basic property tr \ a \ = 2$ab- 
The 't Hooft determinantal interactions are described by the Lagrangian 

£ H = *(det qP L q + det qP R q), (2) 

where the matrices Pl.r = (lT7s)/2 are projectors and the determinant is over flavour indices. 
The coupling constant k is a dimensional parameter ([k] = GeV -5 ) with the large N c asymp- 

totics k ~ 1/N^ S . The coupling G, [G] = GeV~ 2 , counts as G ~ l/iV c and, therefore, the 
Lagrangian dominates over Cu at large -/V c . It differs from the counting G ~ which 
one obtains in the instanton-gas vacuum |15) . 

It is assumed here for simplicity that interactions between quarks can be taken in the long 
wavelength limit where they are effectively local. The 't Hooft-type ansatz (J2J) is a frequently 
used approximation. Even in this essentially simplified form the determinantal interaction has 
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all basic ingredients to describe the dynamical symmetry breaking of the hadronic vacuum and 
explicitly breaks the axial £7a(1) symmetry |3UM31| . The effective mesonic Lagrangian in leading 
order SPA, corresponding to the non-local determinantal interaction, has been found in |32l I33| . 

Anticipating our result, we would like to note that if the hierarchy of multi-quark interactions 
really occurs in nature, the perturbative treatment seems adequate. The NJL interaction alone 
has a stable vacuum state corresponding to spontaneously broken chiral symmetry. But, as we 
shall show, the effective quark theory based on the Lagrangian 1 

£njlh = qii^d^ - m)q + £ NJL + £ H (3) 

has a fatal flaw: if Cu is comparable with £njLj it has no stable ground state. In equation © 
the current quark mass, 777,, stands for the diagonal matrix diag(77\ u , rhd,rh s ), which explicitly 
breaks the global chiral SUl(3) x SUr(3) symmetry of the Lagrangian. 

To bosonize the theory one introduces auxiliary bosonic variables to render fermionic vertices 
bilinear in the quark fields. This procedure requires twice more bosonic degrees of freedom than 
necessary |20j . Redundant variables must be integrated out and this integration is problematic 
as soon as one goes beyond the lowest order stationary phase approximation |34M21| : the lowest 
order result is simply the value of the integrand taken at one definite stationary point [3U]. We 
shall show how to extract systematically the higher order corrections which contribute to the 
effective mesonic Lagrangian, and what to do with infinities contained in these corrections. 

The many-fermion vertices of Lagrangian (J3J) can be linearized by introducing the functional 
unity |U] 



1 



(4) 



/ rjfs a Pp a 5(s a - q\ a q)5(p a - qi^\ a q) 
= J Y[ VsaVp a V<j a V4) a exp |z Jd 4 x [a a (s a - q\ a q) + 4> a (p a - qij 5 X a q)} 
in the vacuum-to-vacuum amplitude 

VqVq exp ( % / d 4 x£ N JLH J • (5) 



We consider the theory of quark fields in four-dimensional Minkowski space. It is assumed that 
the quark fields have colour (N c = 3) and flavour (Nj = 3). The auxiliary bosonic fields, a a , 
and, (j) a , (o = 0, 1, . . . , 8) become the composite scalar and pseudoscalar mesons and the auxiliary 
fields, s a , and, p a , must be integrated out. 

By means of the simple trick (J3J), it is easy to write down the amplitude © as 



/ 8 8 8 8 / f 

VqVq H Vs a JT Vp a JT Va a J] exp (i d 4 xC 
a=0 a=0 a=0 a=0 ^ J 



(6) 



with 



£' = q(i^d^ - a - i^4>)q + s a (a a - rh a ) + p a (j) a + £^JL + £'h> 

^NJL = I [(Sa) 2 + (Paf] , 



£h = T7 I det ( s + + det ( s ~ w)] = -^A abc s a (s b s c - SpbPc) ■ 
54 62 

We assume here that a = cr a A a , and so on for all auxiliary fields a,4>,s,p. 



1 The other multi-quark terms have been neglected here. 
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The totally symmetric constants A a bc are related to the flavour determinant, and equal to 

A a bc = 7^£ijk£mnl(.^a)irn(^b)jn{^c)kl 



2 [2 

-^dabc + Y - (3^a0^60^c0 ~ $a^bc ~ ho^ac ~ ^cO^afe) ■ 



We use the standard definitions for antisymmetric f a bc and symmetric d a bc structure constants 
of U(3) flavour symmetry. One can find, for instance, the following useful relations 

feac-A-bfc + febc-Afac + fefc^-abc = 0, 

deacAbfc + d e b c Af ac + d e f c A abc = V65 e0 A ab f, 
8 [2 8 8 

/ J A-abb = —2y o $a0, A ace A bce = g$ab- 

b=Q c,e=0 

At this stage it is easy to rewrite equation ©, by changing the order of integrations, in a form 
appropriate to accomplish the bosonization, i.e., to calculate the integrals over quark fields and 
integrate out from Z the unphysical part associated with the auxiliary bosonic variables (s a ,p a ) 



Z = J ^Vo a T><^ a VqDq exp (i J d A xC q (q,q,a,( 

X J Yl^SaVpa exp J d A x£ r (a,4>,s,p)^ , (7) 

where 

C q = qii'fdp -a- i~f 5 <t))q, 

C r = s a (a a - rh a ) + p a (f> a + C' N3L + £'h- 

The Fermi fields enter the action bilinearly, thus one can always integrate over them, since one 
deals with a Gaussian integral. One should also shift the scalar fields a a (x) — > cr a (x) + rn a by 
demanding that the vacuum expectation values of the shifted fields vanish ^0|<r a (x)|0) = 0. In 
other words, all tadpole graphs in the end should sum to zero, giving us the gap equation to fix 
the constituent quark masses m a corresponding to the physical vacuum state. 
The functional integrals over s a and p a 



Z[a, eft; A] = N \ ^VsaVpa exp li / d A xC r {a + m, (f), s,] 

J— oo n \ J 



P)) (8) 

are the main subject of our study. We put here A a = m a — m a , and J\f is chosen so that 
Z[0, 0; A] = 1. 

Let us join the auxiliary bosonic variables in one 18-component object Ra = (R a ,Ra) where 
we identify R a = s a and Ra = p a ', o,, a run from to 8 independently. It is clear then, that 
R\ = + Pa- Analogously, we will use ILa = (a a ,(j) a ) for external fields and Aa = (A a ,0). 

Next, consider the sum <&abcRaRbRc- If we require 



$abc = —A a bc, $ abi = — tttAi&c, &abc = $ a bc = ^' 



we find after some algebra 

K 



^abcRaRbRc = C'h 
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with the following important property to be fulfilled 

$abcSbc = 0. (9) 
Now it is easy to see that the functional integral (JSJ can be written in a compact way 

Z[U,A] =M j + \\_VR A exp j d 4 x£ r (n, A; R^j , (10) 



where 



G 2 K 



C r = R A (U A + A a) + -R\ + y ^abcRaRbRc- (11) 
We have arrived at a functional integral with a cubic polynomial in the exponent. 

3 The integration over the auxiliary variables Ra 

Given the cubic structure in the functional integral, one might be tempted to solve it using 
Airy's integral methods. Invoking the existence of a large expansion parameter, such as N c , one 
uses the well known asymptotics of the Airy's function on the real axis. This would require |25| 

C AB = G5 AB + k^abcRc = (12) 



dR A dR B 

which cannot be fulfilled. 

On the contrary, the system of equations based on the first order derivatives 

= GR A + A A + U A + - QabcRbRc = (13) 

is self-consistent and can be solved [2^. Therefore, one can obtain the semi-classical asymptotics 
through the stationary phase method. 

3.1 Solving equation (|13|) 

We need to recall shortly the solutions of equation Up to some order in the external 

(i) 

mesonic fields, n^, we may write them as a polynomial Ra = TV~a 

= Hf + H%U B + H% C U B U C + H A % CD U B U C U D + • • • , 
where i = 1,2,... denote different possible solutions. The coefficients depend on A a and 

(i) 

on the coupling constants G, k, and the higher index coefficients H A are recurrently expressed 
in terms of the lower ones. For instance, we have 

{Hab)' 1 = ~ ( G6 AB + K^ABCH^) , 
ABC ~ 2 ® defH DA H EB H FCi 
and so on. 

Putting this expansion in equation (|13|) one obtains a series of self-consistent equations to 
determine . The first one is 



GHf +A A + | ^ABcHfH^ = 0. 
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One can always find the trivial solution Ha = 0, corresponding to the unbroken vacuum 
A a = 0. There are also non-trivial ones for the scalar component, i.e., 

i?« = (^),0). (14) 

(i) 

The number of possible solutions, i, depends on the symmetry group. The coefficients h a 
are determined by the couplings G, k and the mean field value A a . In accordance with the 
pattern of explicit symmetry breaking the mean field can have only three non-zero components 
at most with indices a = 0, 3, 8. If two of the three indices in A a t, c are from the set {0, 3, 8}, then 
the third one also belongs to this set. Thus A a is the only object which determines the vector 
structure of the solution h a , and therefore h a ^ if a = 0, 3, 8. It means that in general we 
have a system of only three equations to determine = h a ^ X a = diag (hu\ h d \ hi ) 

Gh u + A u + ^-h d h s = 0, 
lb 

Gh d + A d + ^-h s h u = 0, 
lb 

Gh s + A S + ^h u h d = 0. (15) 
lb 

This system is equivalent to a fifth order equation for a one-type variable which can be solved 
numerically. For two particular cases, when m u = rh d = m s and rh u = rh d ^ m s , equations (|15|) 
can be solved analytically The simplest example: m u = m d = m s (or, equivalently, 

hu = h d = hs ) corresponds to SU(3) flavour symmetry. In this case equation (J15|) has two 
solutions 



If AG 2 > kA, they are real and will contribute to the stationary phase trajectory. 

3.2 The lowest order semiclassical asymptotics and instability 

Since the system of equations l(T3]) can be solved, we may replace variables Ra — ► Ra — Ra—T^a 
in the functional integral Q1U|) to obtain the semi-classical asymptotics 

Z[U, A] ~ AA^exp (i J d*xC^ J \\VR A exp Q J ^xC'ab{^)R a Rb^ 

v ®abc J ^xRaRbRcJ (h - 0), (17) 



X 

k=0 



where n is the number of solutions, 71 A , of equation lfT5|) . C" AB has been defined in equation (|T2)) . 
and 



= n$(n A + A a) + ^ (7$) + 1 ^^'4 
= ? fa!?) 2 + 1 + ^) = ^V a + o(n 2 ). 



(0 



Here we used equation (|13|l to eliminate the term proportional to k. Let us stress that C 

(i) . (i) 

depends on k implicitly: k is contained in 7Z A , or more precisely in the coefficients H\ which 
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(i) 

are functions of h a ■ This dependence is singular at k — ► 0. One can see this, for instance, from 
equation (|16() where h u — > oo for small «. This behaviour reflects the fact that we are far from 
the perturbative regime, meaning that the interactions £njl and £h are equally weighted. 

The linear term in the a field is written explicitly. This part of the Lagrangian is responsible 
for the dynamical symmetry breaking in the multi-quark system and taken together with the 
corresponding part from the Gaussian integration over quark fields in equation (J7J) leads us to 
the gap equation. 

At leading order, k = 0, and for the terms linear in the fields contributing to the phase of Z 
we have the estimate for equation (fT7|) 



2 ~ Yl A( S) exp J d * x £-st^J ~ ex P y J d ^ x ~Yl h a )a a + 



where An\ is real and proportional to 2 

Therefore, if one considers the case with m u = rhd = rh s , Z is given by 

Z ~ exp N J d 4 x^ h a )(J a + ■ ■ ■ 1 ~ exp [-i^ f d 4 x(a u + a d + a.) + • • • ^ . (18) 

Let us recall that the quark loop contribution to the gap equation is well known (see, for 
instance, [HS])- Combining this known result with the estimate (|18|). one can obtain the corre- 
sponding effective potential U(m) as a function of the constituent quark mass m 



12G 3N, 
U (m) = m 



16vr 2 



9 

m 



(A 2 -rn 2 ln(l + §)) + <ln(l + ^ 



(19) 



where we consider the case rh = for simplicity and A q denotes the cutoff of quark loop 
integrals. This system has at most a metastable vacuum state (for G/k > 0); for G/k < 0, 
the effective potential does not have extrema in the region m > 0. We must conclude that the 
model considered has a fatal flaw and can be used only in the framework of the perturbative 
approach, which does assume the hierarchy of multi-quark interactions. 

4 Perturbative expansion of Z 

We shall restrict ourselves in this section to the perturbative treatment of the functional integ- 
ral (|lUj) . The loop expansion will be considered in the next section. 

4.1 The perturbative series 

Let us divide the Lagrangian (|11|) in two parts. The free part, Co, is given by 

C Q (R A ) = ^R 2 A + R A (U A + A A ). 
The 't Hooft interaction is considered as a perturbation £j 

Ci(Ra) = j ] <S>abcRaRbRc- 



2 The term constant in the fields of is identical for both critical points. 
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Thus the perturbative representation for the functional integral 1)10(1 can be written as 
Z = N' exp (i j d*x£i(X A )j J l[VR A e * I , ( 2 0) 



where 

S 



X A = ~i 



Since the boson fields appear quadratically in equation (|2()|l. they may be integrated out, 
yielding 

Z = Nexp(iJ tfxd{X A )\ exp f-i J d 4 x^ , (21) 

where fl^ = Ha + Ayj.. The overall factor N = (—2iri/G) 9 N' is unimportant in the following. 

We want to calculate the effective action r c g- = f d 4 x£ e ff> which by definition is the phase 
ofZ 

Z = A(tl A ) exp (*Teff(iU)) , (22) 
and yl(ITi) is a real function. Comparing (|21|) and l|22|) . one gets 

r cff = i In ^ + T - iln (l + e- ir ° ( e ^ d4 ^ - l) e * r «) . (23) 
Here To represents the leading order result for T e ff 

r = ~j d 4 xn A (24) 

while the second logarithm in equation (|23|) is a source of ?7a(1) breaking corrections which arise 
as a series in powers of the functional derivatives operator 



fj = J d 4 xd(X A ) = ^abc J d 4 xX A X B X c . 
To make this statement more explicit let us consider the expansion 



ml 

m=l 



Taking into account the symmetry properties of the coefficients $abc an d our previous re- 
sult ©, we find 



e~ lL 



so that 

oo 

5 = K n S n , 
n=l 

represents the effective action ()23|) as a perturbative series in powers of k. For instance, up to 
and including the second order in k we have 

r cff = i In + T - «<5i - iK 2 f <5 2 - ^5 2 j (25) 
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where 

61 = 31CP ® ABC J d ^ xU ^ U B u c, (26) 

+ m J d 4 x^ + [5(0)] 2 J d 4 ^ . (27) 

The real factor A (II) is always chosen such as to cancel the imaginary part of the effective 
action. To the approximation considered we have, for instance, 

A(ft) = n exp j*(o) (^) 2 J d*x& A } = N (i + 5(0)^ J d 4 xni + •••). 

It contributes to the measure of the functional integral over a a , 4> a . 
4.2 Giving physical meaning to the singularities 

The terms with the 5(0) function require further explanation. The fact that auxiliary fields 
can lead to special problems with infinities is well-known |361 1371 1381 139j . To see in our case 
the origin of the encountered singularities and endow them with physical meaning we shall use 
the language of Feynman diagrams. The graphs contributing to lowest order in n are shown in 
Fig. 1. 





Figure 1. The lowest order graphs contributing to 6±. 

In these diagrams, a line segment (straight or curved) stands for a "propagator" 

A AB (x-y) = -(i/G)S AB S(x-y), (28) 

extracted from the last exponent in equation (|21|). A filled circle at one end of a line segment 
corresponds to the external field, i f d 4 xn^(x), and a vertex joining three line segments is used 
for in&ABC f d x. 

The first diagram represents the 5\ term in equation l|26|l. The contribution of the second 
tadpole diagram is equal to zero. Indeed, the vertex contains the group factor <&abc- The 
contraction of any two indices in this factor by 5ab from the propagator (situation occurring 
for the tadpole graph) reduces it to zero, according to equation @. We thus find that tadpole 
diagrams do not contribute due to the flavour structure of the 't Hooft interaction. 

To next to leading order in k we have the four graphs shown in Fig. 2. 



X 





• + 



e 



Figure 2. The diagrams of the k 2 order contributing to 62- 
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As we just learned, the third diagram does not contribute. The other ones correspond exactly 
to the three terms of equation (|27|). The first tree diagram is finite. The second one- loop diagram 
has a divergent factor 5(0). The last two-loop diagram contributes as 5(0) 2 . These singularities 
were caused by the local structure of the last exponent in equation or, which is the same, 
by the 5(x — y) term in the propagator ()28|) . We believe that if one would start from non-local 
NJL interactions, the singularities could be weaker or would even disappear. 

The factor 5(0) requires a regularization. This is an expected trouble in the NJL model which 
is nonrenormalizable and, as a consequence, the fundamental interactions must be cut off. The 
cutoff is an effective, if crude, implementation of the known short distance behaviour of QCD 
within the model. The problem with 5(0) singularities can also be analysed using the spectral 
representation method to evaluate the integral (jHJ) [25] • For a one-dimensional field theory it 
takes the form 

where L refers to the box size in which the system is put. Therefore the second term does not 
contribute in the limit L — ► oo. The cutoff A cuts the density of Fourier harmonics related 
to the auxiliary bosonic variables, its finite value must be fixed by confronting the model with 
experiment; it may differ from the value of A q in equation (j!9j) . which is associated with the 
integration over the fermionic degrees of freedom. In this interpretation the model is effectively 
finite, including the higher order corrections 3 . 

The diagrams are a very convenient language to understand another feature related to the 
multi-loop contributions: the phase factor corresponding to the diagram can be simply calcu- 
lated. Indeed, it is easy to see that for any diagram the formula 

E = 21 - 3V 

is fulfilled. Here E is the number of external fields, I stands for the number of internal lines, 
and V is the number of vertices. On the other hand, the number of loops, L, is given by 

L = I -V - E + 1. 

This is because the number of loops in a diagram is equal to the number of 5 functions surviving 
after all integrations over coordinates are performed, except for one over-all integration related to 
the effective action. Every internal line contributes one 5 function, but every vertex or external 
field carries an integration over the corresponding coordinate, and thus reduces the number of 
5 functions by one. 

This result shows that the overall phase factor of a given graph i I_E ~ v = i L_1 is entirely 
determined by the number of loops. In particular, diagrams with an even number of loops 
contribute to the effective action (the argument of Z), while the diagrams with odd number of 
loops contribute to the measure (the modulus of Z), see equation (|22j) . 

5 The loop expansion of Z 

The perturbative series, considered in the previous section, can be resummed. This is done 
by introducing a parameter t in the perturbative representation (J20[) of the functional integral 
according to the substitution: 

G, k, Ha, A-a — ► tG, tn, ULa, tAA- 

3 Actually, if one works with the NJL model, one must choose among several known regularizations. Unfortu- 
nately, the dimensional regularization (DR) cannot be used. The gap equation in this case does not have solutions 
and as a result there is no dynamical chiral symmetry breaking. This is why we cannot simply take advantage of 
the well-known result 5(0) = (in DR) to avoid the problem. 



12 



B. Hiller, A. A. Osipov, V. Bernard and A.H. Blin 



Then equation (|25|) to second order in k 2 becomes 



8G 5 



$ABC$AEF / d 4 xnBncn B nir 



t -v-/y - ~( 8 G 2 ) 2 t 2L v ;J 7 32G 3 

The groups in terms of inverse powers of t of T' eS (t) define the new series: first all diagrams 
with no closed loops (tree graphs) will carry the weight 1. The tree graphs yield the result 
of [20]. All the one-loop diagrams (and higher odd numbers of loop diagrams) contribute to the 
imaginary part of the action, as discussed in the previous section; as mentioned before, they will 
be cancelled by the appropriate choice of the real quantity jA(ILa), with A being of order 1. 
They contribute therefore to the measure of the integral over U A and were obtained in |21j . 
see also equations (|3*Tj) . (|3*2*)) below. The set of graphs with 2n loops have the factor t~ 2n and 
contribute to the effective action; in general they include, as a subset, all graphs of fe n+1 th 
order or higher in this coupling constant. The resummed series is equivalent to the well-known 
loop expansion The latter can be written formally as a SPA integral which includes only 
the regular critical point (SPAr), and constitutes the most economic way to obtain the loop 
corrections. A detailed derivation of the equivalence of the resummation of the perturbative 
series, the loop expansion and the SPAr is given in Section 3 of |25| for a one-dimensional 
analogue of the considered functional. 

We proceed therefore to obtain the NLO to the effective action through the SPAr method. 
Consider the effective mesonic action generated by the functional 

+oo 

Z[II,A] -AAexp^i J d 4 x4 =1) ) J l[VR A exp Q J d 4 xC' AB (R^)R A R B ^J 



-co 



00 1 / f \ n 

- ( % ®abc / d 4 xR A R B R c ) (29) 



oo 

X 

n=0 



which in comparison with equation (|17|) has only one critical point, related to the stable 
configuration, which is the solution with i = 1 in equations ()16|) . We shall identify £st _ ^ = C s t, 
and C" AB {T^ i= ^) = C" AB in the following. 

By replacing the continuum of spacetime positions with a discrete lattice of points surrounded 
by separate regions of small spacetime volume fi, the functional integral (|29|) may be reexpressed 
as a Gaussian multiple integral over a finite number of real variables R A (x) for a fixed spacetime 
point x. We think of VR A as the infinite product T>R A — > Y[ x dR A (x), f d 4 x — ► SI Ylx- 

After evaluating the Gaussian integrals one obtains 

• 2 f 



(30) 



Z[U, A] ~ AATT ( Jo exp (iOCO ( 1 + i - ^bcZdefSabcdef . . 
1,1 }\ V 72ttN(N + 2)(N + 4) 

Here 

M = (tr/J^ 1 ) 3 + 6tv£"-hT (C"' 1 ) 2 + 8tr {C"~ l )\ 

contains the complete information on the two-loop term and the dots mean the terms corre- 
sponding to the three-loop contribution and higher. The Jq is the one-loop contribution 



1 /2vr 



- 



K / N 



2 



7 V^7 exp \ 4 ^ 
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where Xj are eigenvalues of the N x N matrix C" AB . In our case N = 18, related to the 
18-component object Ra, equation (|12|). and denotes the inverse matrix of £". The 

totally symmetric symbol Sabcdef generalizes an ordinary Kronecker delta symbol Sab by the 
recurrent relation 

Sabcdef = ^ab^cdef + &ac&bdef H 1- 5af$bcde- 

Up to the given accuracy we have in equation (|3l7|) 

( h e lQCst (1 + iF)) = H ( / e ^ +F ^ 

X X 

= (ll 1 ^ (j]e^ +F ^ = ^/oj e^E.C^+^/n) (32) 

with the quantity Jo containing the one-loop corrections to the measure. It shows that in the 
continuum limit the two- loop correction contributes to the effective Lagrangian as 

C -C 3" 2 W)] 2 M 

LcS ~ Lst + 32N(N + 2)(N + 4)- [M) 

This is our final expression for the effective Lagrangian in the two-loop approximation. 

The field-dependent factor M contains all possible mesonic vertices, including the cr-tadpole 
contribution to the gap equation, contributions to the masses of scalar and pseudoscalar nonets, 
as well as interaction terms. Its dependence on the parameters k, G, m, rh enter through the 
expansion coefficients , equations (|14j) , (|T5|) . 

Let us do some estimates to justify the result. For this purpose let us simplify the integ- 
ral (|29|). After neglecting the symmetry group and discretizing the spacetime it takes the form 

Z[U, A] ~ / dRx exp (c st + ^ t R 2 x + } . (34) 

To justify the stationary phase approximation for the integral (|3"1|) we assume that 

nC st > 1. (35) 

The dominating role of the Gaussian integral is reflected in the fact that essential values for R x 
in the integral have the order R% ~ l/(0£g t ). For the cubic term it follows then 



where we have used that in the model considered here, C"l ~ k, C" t ~ G. If the parameters of 
the model can be chosen in such a way that the inequality £ <S 1 is fulfilled, the cubic power 
of R x yields terms that go to zero relative to the Gaussian term as £ — > 0, and the stationary 
phase approximation will be justified. Note that may be written as an ultraviolet divergent 
integral regularized by introducing a cutoff A 



Q- 1 = S\0) 



A / 2 d*kE = (AY 
^ A/2 (2vr) 4 ~ \27t) 



Therefore, the inequality restricts the value of A from above. 
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Meanwhile, it is interesting to see that the inequality £ <C 1 is an exact equivalent of (13511 . 
Indeed, in the essential region, i.e., around a sharp minimum, one has C' st ~ R x C s t, C" t ~ R%.C s t, 
and so on, thus 



c 



m\2 
st) 



{RlC 



st.J 



n(£«)3 mi^st) 3 oc st ' 



Therefore, the asymptotical series (|30|) with the ultraviolet cutoff imposed in the contin- 
uum limit is sensible. One deals here, actually, with a series in powers of the dimensionless 
parameter £. The expansion is formally justified for 

5.1 The NLO effective potential 

As an example let us obtain the NLO corrections due to (|33|) to the effective potential for the 
SU(3) limit and with the current quark masses set to zero. Here only the regular critical point 
was taken into account. For the case of the full SPA, one has to consider the same expression also 
for the singular point. The term linear in the scalar fields coming from the two-loop correction 
is calculated with the help of "Mathematica" j3Uj 

2 Yl Ci > a(Ta = + a d+ <T t ) 



(121 - 1552a; 2 + 7832^ - 17980a;? + 14467wf + 3592o^°) 
X (l-4u; t ?)4(l-a; l ?)4(l + 2a; i ) 

and it must be added to the leading order term equation ((T5|) . Here the coefficient 



A 
2^ 



3k 2 



32N(N + 2)(N + 



and tOi = L0 U) i = L0d,i = ^s,i = with the index i denoting the two possible solutions for hi to 
the SPA equations in the SU(3) case, equation (|16|) . One must integrate q with respect to the 
quark mass to obtain the correction to the effective potential, since 

E8G 2 
Ci — (1 + 2iOi)dwi, 



dV„ 



where for the last equality the equations (|15j) were used to obtain 

16G 2 



dm(i) = 
Here m = m u 

C^lcorr — 



f ^ ) dhi 



m d 



48 



-(1 + 2uji)(hji. 

m s and m(i) = m for both solutions. With X{ = cof one gets 
121 - 1552a;; + 7832x 2 - 17980a; 3 + 14467xf + 3592xf ) 



(l-4x i ) 4 (l-x i ) 4 



dxi. 



As a result we obtain 



V a 



IS 



55 - 583xi + 2023a; 2 - 2037k 5 - 898xf 



G 3 'Y 2(l-x l ) 3 (l-4x i ) 3 
Therefore the complete expression for the SU(3) effective potential is 
12Gm 3N C T o / . n „ / A 2 



U(m) 



16vr 2 



m 2 A 2 - m 2 In 1 + -J 



A^ln(l + ^ 



771 
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There will be poles at Xi = {j, 1}. In terms of m there can be in principle up to 4 non-degenerate 
poles altogether, two for each critical point. They are due to the zero eigenvalues of C", which 
appear for certain values of the coupling constants |21| in subsets of the 18-field configurations. 
In the neighborhood of these points the SPA method is unreliable. This puts constraints on 
the coupling constants, which must be chosen such that the poles appear far from the region of 
physical interest (the quark mass value at the local minimum of the effective potential). One 
sees however that the solution /12, related with the singular critical point, has always a pole at 
m = 0, since for this case X2 = 1. It is one more indication that the contribution of the singular 
point must be excluded from the analysis. Only the perturbative regime is sensible. In this case 
the effective potential reads 



3Ghj + 



3Nc 
16vr 2 



m 



A 



m 2 ln f 1 + ^| 



1+ Af) 



+ V^ COTT (m) 



where the index 1 labels the contribution of the regular point. The solution h\ will lead to poles 
at m = 12 p and m = 32 . The parameters G, k can be chosen such that the pole positions 
are far away from the constituent quark mass solution of the gap equation. For instance for 
« = -1300 GeV" 5 , A = 1.8 GeV, A„ = 0.82 GeV the two poles due to hi 



G = 13.5 GeV -2 , k 

appear at m ~ 1.7 GeV and m ~ 4.5 GeV. For G = 2.15 GeV" 2 , k = -53 GeV" 5 , A = 1.9 
GeV, A q = 1.64 GeV one has the two poles at m ~ 1. GeV amd m ~ 2.8 GeV. Both sets 
of parameters lead to good meson mass spectra and weak decay constants in the realistic case 
with SU(2) x C/a(1) symmetry |41) . In conclusion, from the point of view of the complete SPA 
(inclusion of regular and singular points), there is no improvement, because of the bad pole at 
zero. From the point of view of the perturbative (regular) solution, the poles are far away on the 
positive m axis and the present approach can be used to study systematically NLO corrections. 
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